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Disclaimer 
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expressed  or  implied,  of  the  Advanced  Research  Projects 
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Summary  of  Research 

The  primary  purpose  of  this  research  was  to  develop  the 

phenomenon  of  extended  x-ray  absorption  fine  structure  (EXAFS)  as  a 

technique  for  investigating  the  structure  of  non-crystalline  materials, 

particularly  amorphous  semiconductors.  A  theory  had  been  developed 

(Sayers,  Lytle,  Stem,  "Advances  in  X-ray  Analysis",  Vol.  13,  p.  2A8 

(Plenum,  New  York,  1970)  which  treated  EXAFS  as  a  scattering  of  the 

ejected  photoelectron  by  the  near  neighbors  of  the  absorbing  atom. 

Soon  after  the  research  began  a  major  breakthrough  was  achieved  in 

which  the  theoretical  EXAFS  expression  and  xperimental  data  could  be 

formally  inverted  into  a  radial  structure  function  using  the  Fourier 

Integral  theorem The  technique  was  further  developed  and  applied 

(2) 

to  an  investigation  of  the  structures  of  Ge,  GeO^.  and  GeSe  to  make 

the  important  point  that  the  environment  of  all  the  different  atoms  in 

a  multicomponent  material  could  be  determined.  This  work  was  recognized 

as  a  significant  contribution  to  x-ray  physics  when  Dale  Sayers  was 

awarded  the  1972  Sidhu  award  by  the  Pittsburgh  Diffraction  Conference. 

Dr.  Sayers'  Fh.D.  thesis,  which  described  this  wurk  in  detail  was  sub- 

(3) 

mitted  as  a  special  technical  report. 

The  surface-active  character  of  evaporated  amorphous  germaniun  (E  Ge) 

(4) 

was  emphasized  in  a  calorimetric,  diffraction  and  microscopic  evaluation 
which  showed  that  3  wt.Z  could  be  adsorbed  when  E  Ge  was  exposed  to 
air  and  could  not  all  (50%)  be  removed  unless  the  sample  was  heated 
above  the  re  crystallization  temperature.  During  recrystallization  electron 
diffraction  patterns  showed  transient  (222)  reflections  (normally  forbidden) 
which  were  interpreted  as  evidence  of  wurtzite-type  bonding  in  the 
original  material. 

The  EXAFS  work  on  amorphous  semiconductors  was  continued  on 
GeSe,  GeSe2>  As^e^.  As,^.  and  As2Te^^  and  reported  at  the  1973 
Conference  on  Amorphous  Semiconductors  at  Garmisch-Fartenkirchen. 

This  work  developed  a  normalization  technique  which  allowed  determination 
of  absolute  coordination  numbers  as  well  as  identifying  the  kinds  of 
atoms  in  the  first  coordination  sphere.  For  example,  in  amorphous 


As  Se  we  found  the  As  atom  to  have  3  nearest  neighbors  (2Se  and  1  As) 

2  3 

while  Se  had  2  nearest  neighbors  (1  As  and  1  Sc). 

The  technique  was  also  used  to  investigate  impurity  charge 
shielding^  and  determine  atomic  arrangement  around  the  catalytic 
atoms  in  a  supported  catalyst^.  The  technique  is  finding  great 
Interest  among  catalyst  chemists  at  a  time  when  progress  in  this  impor¬ 
tant  energy— related  field  is  in  the  national  interest. 

Accompanying  this  final  report  are  three  papers  <whleh  will  he 
submitted  *£or  publication)  describing  the  EXAFS  experimental^  ^  and 
data  analysis  techniques^9’ These  reports  describe  in  detail  our 
practice  of  the  technique. 

For  the  future  of  the  technique,  a  proposal  was  made  to  NSF 
and  a  grant  awarded  to  the  University  of  Washington  (Stem,  Sayers, 
Lytle-consultant)  to  develop  an  EXAFS  experiment  at  the  Stanford  Synchro¬ 
tron  Radiation  Project.  The  x-ray  intensity  which  can  be  obtained  at 
this  source  is  at  least  lfi5  greater  than  that  of  ordinary  x-ray  tubes 
and  will  allow  data  of  higher  resolution  and  better  precision  in  a  much 
shorter  time. 

Conclusions 

In  competition  with  x-ray,  electron,  or  neutron  scattering 
the  EXAFS  technique  has  the  important  advantage  that  the  origin  of  the 
radial  structure  function  (two  body  correlation  function)  is  fixed  at 
the  site  of  the  species-determined  x-ray  absorption  event  while  all 
other  techniques  average  over  all  the  different  kinds  of  atoms  in  the 
material.  The  results  obtained  on  amorphous  semiconductors  are 
important  and  unique  and  will  allow  theorists  to  construct  more  accurate 
models  of  these  materials.  The  catalyst  work  is  only  the  top  »f  an 
iceberg.  This  was  the  first  example  of  a  technique  in  whidi  tut! 
environment  of  catalytic  atoms  in  general  could  be  determined. 

The  EXAFS  technique  has  been  shown  to  be  a  powerful  tool  for 
studying  the  arrangement  of  atoms  in  multi-component,  non- cry s tall ine 
materials. 
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*  * 

Farrel  W.  Lytle  and  Dale  E.  Sayers 

The  Boeing  Company,  Seattle,  Washington  98124 

Abstract 

Techniques  of  x-ray  spectroscopy  for  obtaining  extended  x-ray 
absorption  fine  structure  (EXAFS)  are  discussed.  The  essential 
instrumentation  is  an  automated  commercial  x-ray  diffractometer 
with  associated  power  supply  and  counting  electronics  modified  to 
include  an  improved  x-ray  monochromator  alignment  fixture,  a  cryostat 
for  containing  the  absorber  sample  and  a  card  punch  for  data  readout. 
The  preparation  of  suitable  monochromator  crystals  is  discussed. 

A  preset  count  mode  was  used  to  make  the  most  efficient  i»e  Oi 
experimental  time  and  the  technique  of  multi-pass  averaging  was  used 
to  obtain  good  counting  statistics;  i.e.  0.1%.  The  precision  and 
accuracy  of  the  wavelength  /energy  scale  was  determined  by  replicate 
.analysis  and  comparison  to  emission  lines  from  the  x-ray  source  to 
be  leV.  When  these  lines  were  intense  the  nonlinear  response 

I 

of  the  x-ray  detecting  system  required  a  coincidence  correction 
to  be  evaluated  and  applied  to  the  measured  counting  rates.  EXAFS 
data  from  Cu  is  shown  as  an  example. 


*  Research  supported  in  part  by  the  Advanced  Research  Projefcts 
Agency  and  monitored  by  the  U.  S.  Army  Research  Office, 
Contract  No.  DAH  C04  71  C  0010. 
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Introduction 


This,  and  the  following  paper,  describe  the  experimental  details 
and  data  processing  techniques  which  we  have  developed  for  the  prac¬ 
tice  of  x-ray  absorption  spectroscopy  for  determination  of  the  struc¬ 
ture  of  non-crystalline  solids  (For  a  genera^  treatment  of  x-ray 
spectroscopy  techniques,  see  Sandstrom^  .)  In  our  technique  the 
extended  x-ray  absorption  fine  structure  (EXAFS)  is  Fourier  analyzed 
to  produce  a  radial  distribution  function  centered  on  the  absorbing 
atom^.  The  particular  utility  of  the  technique  is  the  ability  to 

xplore  the  average  environment  of  each  atomic  species  in  multicom- 

(3) 

ponent  materials  such  as  amorphous  semiconductors  or  the  highly 

(4) 

dispersed  materials  used  as  catalysts  .  The  basic  experimental 
determinant  of  the  technique  is  a  high  resolution,  precise  measure¬ 
ment  of  the  x-ray  absorption  coefficient  in  the  neighborhood  of 
characteristic  absorption  edges.  This  paper  describes  the  use  of 

a  commercial  diffractometer  (with  fairly  simple  modifications)  to 

* 

produce  an  x-ray  spectrometer  of  moderate  resolution;  the  diffrac¬ 
tion  crystals  and  preparation  of  the  diffractive  surfaces;  x-ray 
detectors  and  counting  techniques;  and  a  cryostat  and  sample  trans- 
port  mechanism. 

Summary  of  theory 

During  an  atomic  x-ray  absorption  event,  an  x-ray  with  energy 
ho  disappears  and  a  photoelectron  of  kinetic  energy  E  is  emitted 
from  the  atom  conserving  energy  in  the  process 

E  -  hu  -  Er  ll) 


where  is  the  initial  binding  energy  of  the  electron.  The  absorp¬ 
tion  follows  the  law 

ux  ■  In  I  /I  or  In  T/T  (2) 

.  ~  o  o 

where  is  the  linear  absorption  coefficient,  x  is  the  thickness, 

I  is  the  incident  x-ray  intensity,  1  is  the  transmitted  intensity 

a 

and  T  and  Tq  are  the  times  to  preset  count  with  the  absorber  in 

and  out  of  the  beam.  On  the  high  energy  side  of  the  absorption 

/ 

edge  px  shows  fluctuations  (EXAFS)  which  are  related  to  the  atomic 

arrangement  immediately  surrounding  the  absorbing  atom.  We  have 

shown  that  the  EXAFS,  X(k)  can  be  accounted' for  by  a  scattering 
(5) 


model 


N 


(3) 


X(k)  -  hY.^2  E(r  )g  (2k)  exp(-2o  V)  Sin[2kr  *  +  2T)  (k)] 
j  J  J 

where  k  is  the  electron  wave  vector,  is  the  number  of  atoms  in 

the  1th  coordination  sphere,  rA  is  the  radial  distance  to  the  jth 

j 

atoms,  E(r  )  is  the  scattering  potential  encountered  by  the  electrons, 

«f 

g  (2k)  is  the  backward  scattering  amplitude,  the  exponential  is  a 

J 

Debye-Waller  type  term  where  (J,  is  the  rms  fluctuations  of  the  atom 

•J  / 

about  r^  and  r|(k)  is  a  phase  shift.  The  form  of  this  equation  makes 

»  EXAFS  proportional  to  the  number  of  atoms  surrounding  the  absorbing 

2 

atom  and  inversely  proportional  to  r^  .  A  decrease  in  temperature 
has  its  usual  effect  of  sharpening  the  EXAFS.  Each  coordination 
sphere  contributes  a  sine-like  term  of  period  2kr ^ . 

Conceptually  EXAFS  may  be  visualized  as  a  kind  of  electron 
diffraction  where  the  source  of  the  electrons  is  within  the  material 
originating  from  the  atoms  of  that  particular  atomic  species  involved 
in  the  absorptiou  event.  The  mechanism  of  the  absorption  fluctua¬ 
tions  is  as  follows:  during  the  dipole  transition  the  electron 
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oscillates  between  the  initial  and  final  states.  Since  the  final 
state  in  this  case  is  located  away  from  the  absorbing  atom  the  elec¬ 
tron  encounters  and  is  scattered  by  .he  surrounding  potential  array. 

The  overlap  of  the  scattered  and  initial  wave  functions  in  the  dipole 
transition  matrix  produces  the  EXAFS  fluctuations. 

X-ray  Spectrometer 

Our  basic  x-ray  spectrometer  (as  shown  in  Fig.  1)  is  a  hori¬ 
zontal  diffractometer  with  attached  tube  stand.  The  only  modifi¬ 
cation  is  an  improved  crystal  (monochromator)  support  for  optimum 
alignment  of  the  diffracting  crystal  (Fig.  lb).  This  spectrometer 
is  used  with  conventional  x-ray  diffraction  tubes  (Ag  or  Mo  target) 
chosen  to  provide  a  high  intensity  continuum  in  the  region  of  the 
absorption  edge  *■ o  be  measured.  The  problem  of  intense  spectral 
lines  which  sometimes  occur  is  discussed  in  the  section  on  counting 
techniques.  We  use  the  spectrometer  with  a  goniometer  radius  of 
18.75  cm  and  an  x-ray  tube  takeoff  angle  of  approximately  3  (adjust¬ 
ed  for  maximum  diffracted  intensity)  .  The  first  slit  is  used  to 
limit  the  angular  divergence  of  the  x-ray  beam  (tube  spectrometer 

l 

mode)  and  the  Bragg-Brentano  parafocusing  conditions  and  slit  posi¬ 
tioning  require  the  focusing  (exit)  slit  to  be  approximately  thr 
same  size.  The  vertical  divergence  is  limited  to  ±  2°  by  a  Soller 
slit  located  before  the  exit  slit.  Thus,  the  resolution  function 
depends  primarily  on  the  size  of  the  entrance  slit.  (See  Klug  and 
Alexander for  a  complete  discussion  oi  this  type  of  spectrometer.) 
Our  usual  method  of  operation  employs  0.05  mm  slits  (.025  divergence) 
in  the  range  from  15°  to  45°  20  and  0.1  mm  (.05°  divergence)  slits 
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for  angles  greater  than  45°  with  overlap  as  experimentally  required. 

Hence,  the  resolving  power  defined  as  0/A0  is  approximately  1500. 

The  diffracted  intensity  from  the  Continuum  operating  at  full  recom- 

3 

mended  tube  KVA  depends  on  wavelength  but  is  typically  1-10x10  cps. 

For  the  case  of  0.05  mm  slits  at  45°  20  and  a  LiF  crystal, 
the  energy  band  width  received  at  the  exit  slit  would  be  ~3.5eV. 

Since  the  intensity  distribution  of  the  radiation  filling  the  exit 

i 

slit  has  the  usual  diffraction  profile,  the  Rayleigh  resolution 
criterion  suggests  that  the  spectrometer  should  be  advanced  in  angu¬ 
lar  increments  of  *s  the  angular  width  of  the  exit  slit. 

The  mode  of  operation  is  as  follows:  For  each  spectrometer 
position, T  and  T  are  measured  (preset  count)  and  stored  on  IBM 
cards,  the  spectrometer  is  advanced  to  the  next  20  position  and  the 
sequence  is  repeated.  A  separate  scaler  is  used  to  generate  a 
running  number  for  each  subsequent  pair  (T,  Tq)  .  Knowledge  of  the 
start  position  and  the  20  stepping  increment  allows  calculation  of 
the  x-ray  wavelength  for  a  ay  data  pair.  The  usual  experiment 
involves  500  -  1000  data  pairs.  The  spectrometer  stepping  motor 
*  and  absorber  changer  are  activated  and  synchronized  by  the  x-ray 

scaler  print-out  command.  • 

The  flanged  tube  stand  attached  to  the  spectrometer  allows  the 
x-ray  tube  to  be  transtated,  rotated,  and  inclined  with  respect  to 
the  spectrometer  circle.  The  alignment  procedure  consists  of 
locating  (by  means  of  these  adjustments)  the  most  intense  spot  on 
the  x-ray  tube  target  so  that  it  shines  through  the  zeroed  entrance 
and  exit  slita  as  measured  by  a  protected  x-ray  detector.  The 
crystal  monochromator  is  then  inserted  and  a  suitable  characteristic 
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line  chosen  for  final  crystal  adjustment.  At  the  calculated  20 
the  translation  and  rotation  (0)  adjustments  of  the  crystal  holder 
are  used  to  obtain  peak  diffracted  intensity  of  the  standard  line 

m 

which  locates  the  diffracting  surface  of  the  crystal  at  the  center 
of  the  spectrometer.  Further  refinement  of  the  alignment  is  not 
necessary  if  the  procedure  described  in  the  section  on  precision 
and  accuracy  of  the  energy  scale  is  followed. 


Cryostat  and  Absorber 


'  *  (7) 

'•EXAFS  has  been  shown  to  be  sensitive  to  temperature  .  The 

degree  of  the  effect  is  illustrated  in  Fig.  2b  and  2c .which  were 

obtained  at  77°K  and  295°K,  respectively.  The  Debye-Waller  type 

• 

exponential  occurring  in  the  theoretical  expression  for  EXAFS  can 

(8) 

be  related  to  the  usual  Debye  expression v  . 

o,2  -  d  A,  (T/en)  ■  C4) 

where  0  is  the  Debye  characterisitic  temperature,  c.  is  the  rms 
D  J 

amplitude  of  vibration  of  the  atoms  in  the  jth  coordination  sphere, 

K  is  Boltzmann  constant,  m  is  the  mass  of  the  jth  atoms,  di  is 
*  B 

Planck's  c  cnstant/2ir,  ^(T/0D)  is  the  Debye  function,  and  T  is  the 
absolute  temperature.  Thus,  for  materials  with  a  low  value  of  0^, 
the  temperature  effect  will  be  largest,  however,  all  materials  will 
show  an  EXAFS  intensity  increase  (particularly,  far  from  the  absorp¬ 
tion  edge).  For  this  reason  data  is  obtained  routinely  at  77  K 
(extra  sharpness  obtained  at  4  K  isn't  worth  the  effort).  It  is 

possible  to  use  the  temperature  effect  to  directly  evaluate  the  shell- 

(8)  • 

to-shell  coupling  coefficients 
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PHOTOELECTRON  KINETIC  ENERGY,  eV 


Figure  2.  Plots  of  In  T/_  vs  the  energy,  E,  of  the  ejected 
photoelectron  for  -2.5pm  thick  Cu  foil. 

Figure  2a.  One  pass  at  77°K,  10^  preset  count.  The  features 

coincident  with  the  W  emission  lines  are  due  to  the 
nonlinearity  of  the  x-ray  detection  system.  This 
effect  has  been  removed  on  the  other  graphs  using 
a  correction  technique  described  in  the  text. 

Figure  2b.  Ten  pas1  as  at  77°K,  10^  preset  count  each  pass. 

Figure  2c.  Ten  passes  at  295°K,  105  preset  count  each  pass. 

Note  the  smearing  of  fine  structure  particularly 
at  high  energy. 
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The  glass  cryostat  Illustrated  In  Fig.  1  contains  5  liters 
of  liquid  nitrogen  which  will  last  for  3  days.  The  design  of  the 
sample  moving  mechanism  (Fig.  1M  accomplishes  Its  purpose  without 


movable  vacuum  seals  or  introduction  of  a  heat  leak  into  the  sample 

r*  6 


area.  If  a  stainless  steel  bellows  Is  used  the  fatigue  life  Is  10 
cycles  or  1-2  years  of  operation.  The  x-ray  transparent  windows  are 
1  all  Mylar  expoxied  to  oval  openings  In  the  outer  shell. 

The  absorber  samples  are  prepared  in  a  variety  of  ways;  malleable 
metals  are  rolled  (2-5Mm) ,  some  materials  may  be  evaporated  onto 
Mylar  or  thin  A1  foil,  soluble  materials  may  be  dissolved  and  then 
absorbed  and  dried  on  filter  paper,  many  materials  may  be  ground  as 
fine  as  possible,  mixed,  with  a  suitable  plastic  glue  (Duco)  and 
solvent  (acetone)  then  cast  on  a  smooth  substrate.  When  dry,  the 
casts  are  sandwiched  between  H  mil  A1  for  support  and  thermal  contact 
and  attached  to  the  sample  holder.  The  optimum  absorber  thickness 
considering  contrast,  measurement  time,  and  primarily  the  coincidence 
counting  error  problem  is  attained  when  I/I.~  1/3  on  the  high  absorp¬ 
tion  side  of  the  edge.  The  experiments  when  only  a  very  small  amount 
'  of  the  material  of  interest  is  present  must  use  thicknesses  such  that 

I/I  -  1/10.  For  example  with  this  thickness  a  satisfactory  Lm 
'  O 

e  n  ’X  uf  X  Pt  on  A1„0„  was  made  in  15  passes.  - 
absorption  pattern  of  0.3  wt  S  rt  on 


Mono rhromator  Crystals 


The  usual  discussion  of  monochromators  for  x-ray  spectrometers 
emphasises  resolution  with  narrow  crystal  rocking  curves  and  multi¬ 
crystal  spectrometers.  Our  requirement  stresses  high  x-ray  intensity 
at  moderate  resolution  for  good  statistical  accuracy  of  measured  EXAFS . 


M* 

- - 
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Although  better  resolution  would  probably  show  more  detail  there 
is  an  inherent  width  in  the  EXAFS  due  to  temperature  smearing  and 
uncertainty  broadening ^  and  a  spectrometer  with  better  resolution 
than  the  natural  width  is  wasted  effort. 

The  x-ray  intensity  is  limited  by  the  inherent  brilliance  of 
the  x-ray  tube  and  using  flat-crystal  optics  the  count  rate  is 
limited  by  the  efficiency  of.  the  diffracting  crystal  and  the  width 

t 

of  the  energy  band  through  the  exit  sHt.  The  diverging  x-rav 

beam  has  some  distribution  in  space  and  energy;  i.e.  photons/sec/ 

stearadian/eV.  The  monochromator  selects  with  some  diffraction 

efficiency  a  band  of  radiation  equal  to  the  angular  width  of  its 

rocking  curve  and  directs  it  toward  the  exit  slit.  The  optimum 

% 

intensity  will  occur  when  the  divergence  of  the  incident  beam  as 
limited  by  the  first  slit  and  the  crystal  rocking  curve  are  appro¬ 
ximately  equal. 

The  efficiency  of  diffraction  and  the  rocking  curve  width  of 
diffraction  crystals  may  be  modified  by  appropriate  surface  treat¬ 
ment.  LiF  is  a  particularly  workable  crystal  in  this  respect.  In 
♦Fig.  3  the  integrated  reflection  coefficient,  R,  for  non-polarized 
radiation  which  is  the  area  under  the  crystal  rocking  curve,  has 
been  calculated  as  a  function  of  wavelength  for  the  two  extremes  of 
crystal  perfection,  a  perfect  crystal  and  an  ideal  mosaic,  crystal^  ^ 
The  measured  values  of  selected  cleaved  and  treated  crystals  are 
also  shown.  Illustrated  in  Fig.  A  a  and  b  are  the  full  width  at 
half  maximum  (fwhm)  of  the  rocking  curve  and  the  diffraction  effi- 

r 

ciencv  (l  ,  ,1  ,  ,  )  at  the  peak  of  the  rocking  curve, 

ciency  Udlff racted/  incident'  v 

The  measurements  were  made  on  a  double  crystal  spectrometer  in  the 
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DEGREES  DIFFRACTION  EFFICIENCY, 


Figure  3,  Integrated  reflection  coefficient  vs  wavelength  for 
LiF  as  cleaved,  roughened,  and  chemical  polished. 

Hie  measurements  were  made  on  a  double-crystal 
spectrometer  correcting  for  the  partial  polarization 
of  the  incident  x-ray  beam. 


Figure  4b.  The  intensity  ratio,  diffracted/incident,  at  the 

peak  of  the  diffraction  curve  for  LiF  as  a  function 
of  wavelength  and  surface  treatment. 


>  wavelength  and  surface  treatment. 


(1,  -1)  mode  using  the  first  crystal  to  select  a  monochromatic 
beam  of  small  angular  divergence  (<  .01°).  The  test  crystal  was 
then  "rocked"  (rotated)  through  its  Bragg  angle  while  the  dif¬ 
fracted  photons  were  measured.  This  is  a  direct  measure of 
R  since 

*  -  fi!_  -  f  lliL.  d0  (5) 

Iin.  J  'in. 

where  £  is  the  total  number  of  photons  counted  under  the  rocking 
curve,  w  is  the  angular  rate  of  crystal  rotation  (in  radians)  and 
I^n  and  refer  to  the  incident  and  diffracted  radiation. 

The  measurements  have  been  corrected  for  the  partial  polarization 
of  the  incident  test  beam. 

These  results  show  that  the  total  radiation  diffracted  by  LiF 
may  be  increased  by  factors  of  2-5  by  suitable  surface  treatment  to 
reduce  the  extinction  effects  which  occur  during  diffractions 
The  treatment  consists  of  vigorous  sarding  on  rough  paper  to  drive 
dislocations  into  the  crystal  followed  by  successively  finer  paper 
to  600  grit  to  provide  a  smooth  surface.  This  was  the  treatment 

*  of  the  abraded  crystal  in  Figs.  3  and  4.  The  damaged  surface 

(12) 

layer  was  then  removed  by  an  etching  procedure  .  After  1  minute 
in  cone.  HF  a  chemical  polish  consisting  of  2%  NH^OH  in  at  26°C 
with  vigorous  agitation  was  used  to  remove  the  surface  at  about 
lpm/mln.  The  crystal  was  repeatedly  checked  until  its  rocking  curve 
width  narrowed  to  that  of  the  slits.  Note  that  for  the  etched  crystal 
the  diffraction  efficiency  and  R  were  considerably  enhanced  over 
the  cleaved  crystal.  We  have  used  similar  techniques  with  Si  and 
quartz  crystals  with  some  success. 


Counting  Technique 


The  precision  to  which  the  EXAFS  function  may  be  determined 
should  be  limited  by  the  usual  statistics,  og  where 

standard  deviation  and  N  -  number  of  photons  counted  for  each 
data  points  This  is  usually  expressed  as  percent  counting  error, 
pee  -  (og/N)  x  100%;  e.g.  when  N  -  106,  pee  -  0.1%.  It  would 
appear  that  any  level  of  precision  could  be  obtained  by  making  N 
sufficiently  large;  but  in  practice  the  inherent  stability  (elec¬ 
trical  and  mechanical)  cf  the  experiment  will  limit  the  precision. 

It  is  important  to  determine  this  value  for  a  particular  experiment 
otherwise  time  will  be  wasted  by  operating  at  a  larger  N  than 
warranted. 

When  the  inherent  experimental  precision  has  been  evaluated 
(as  described  below)  the  experiment  is  set  to  accumulate  the 
appropriate  number  of  counts  for  each  data  point.  If  better  pre¬ 
cision  is  desired  then  multi-pass  averaging  is  used  to  accomplish 

it. 

The  inherent  precision  of  our  experiment  was  determined  as 
'  follows :  For  a  rang,  of  N  from  10-107  a  subset  of  data  was  accumu¬ 
lated  for  each  value  of  N  and  ss  estimated  using_  standard  statistics 

os  -  [Sf  [K/I,)  -  (  VP  ]  2  ]  '  (6>  ’ 

where  Tq  and  T  refer  to  the  time  to  accumulate  N  counts,  the  barred 
value  is  the  mean  of  the  subset,  and  n  is  the  number  of  trials  in 
each  subset.  The  result  of  this  analysis  is  plotted  in  Fig.  5  where 
the  difference  between  the  measured  and  predicted  pee;  i.e. 

•  5 

[o8 (measured)  -  N5]  100/N  is  plotted  vs  log  N.  Up  to  N  = 
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!  Ui 

g  LOG  ACCUMULATED  COUNTS,  N 


Figure  5.  The  difference  between  the  predicted  and  measured 
percent  counting  error  vs  total  counts. 


lo/|  AND  CORRECTION  FACTOR 


Figure  6.  Counting  rate  vs  absorption  ratio,  I Q/I,  for  the 
detecting  systems  a,  b,  c  described  in  Table  2. 

In  the  region  of  linear  response  the  ratio  is  constant, 
then,  since  I  >  I,  I  decreases  relative  to  I  and  the 
absorption  ra?io  decreases.  Graph  d  is  a  plot  of  the 
correction  factor  (Rhort’sC14)  method)  for  nystem  b 
which  can  bo  multiplied  times  the  niconured  count  rate 
to  correct  for  nonlinearity. 
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the  data  shows  a  ra.idoro  scatter  reflecting  the  accuracy  with  which 

a  was  measured  but  at  N  -  105  there  is  a  t^arp  break  in  the  experi- 
s 

mental  curve  indicating  that  a  constant  level  of  precision  has  been 
attained  no  matter  how  much  larger  N  was  made.  If  th.e  experiment 
had  been  completely  statistical' /  limited  the  pee  data  points  would 
cluster  around  the  zero  axis  for  all  N.  It  is  clear  from  this 
experiment  that  accumulating  more  than  10  counts  (0.3%  pee)  is  a 
waste  of  time.  The  reason  for  this  fairly  high  pee  was  probably 
due  to  slight  mechanical  instability  in  the  crystat  support  which 
allowed  the  absorber  to  move  in  the  x-ray  beam  and  reflect  irregu¬ 
larities  in  thickness  in  the  measured  signal.  The  improved  pre¬ 
cision  apparent  in  the  EXAFS  of  Fig.  2b  and  2c  was  attained  by 

averaging  ten  separate  runs  at  10  counts  for  each  data  point  for 

0 

a  final  precision  near  0.1%. 


Precision  and  Accuracy  of  the  Energy  Scale 


The  kinetic  energy  of  the  ejected  photoelectrons,  E,  must  be 
established  accurately  in  order  to  evaluate  the  natural  EXAFS 

t 

variable  k  where 

k  «=  2m^  (E  -Aw)^  **  (7) 

where  Aw  is  the  difference  between  the  work  function  and  the  aver¬ 
age  potential  seen  by  the  scattered  electron.  For  every  experiment 
characteristic  and/or  impurity  lines  from  the  x-ray  tube  occur  and 
can  be  used  as  standard  reference  points  to  establish  tue  energy 
scale  at  the  accuracy  to  which  they  are  known  .  As  an  example, 
the  lines  indicated  in  Fig.  2  which  occurred  in  the  region  of  the 
experiment,  are  tabulated  in  Table  1.  Many  elements  are  present  as 


Table  1.  Precision  of  Energy  Scale 


Characteristic 

Line 

Wavelength, 

1.2047917 

Measured  20 
Average  of  10 runs 

Calculated  2d 

W  L  P9 

34.83813 

4.02459 

W  L  &5 

1.2154513 

35.16212 

4.02396 

W  L  P2 

1.2/46013 

36.02611 

4.02480 

W  L  &3 

t 

1.2626915 

36.56411 

■ 

4.02523 

w  l 

1.23180919 

37.14212 

4.02478 

W  L  P6 

* 

1.2898917 

37.38712 

4.02455 

, 

- 

w  l 

1.3016215 

37.73312 

** 

4.02524 

Cu  K  ^ 

1.39221819 

40.47312 

4.02496 

w  l  n 

1.4211013 

41.35012 

4.02502 

• 

Avg. 4.024713 
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impurities  on  the  x-ray  tube  anode;  e.g.  W(always  present-eval- 
porated  from  the  filament),  Cu  (always  present-evaporat-.d  from 
the  anode  block),  Pb,  Ni,  Ag,  and  Zn  are  usually  present,  plus 
the  characteristic  lines  from  the  primary  anode  element  in 
multiple  orders  of  diffraction.  The  presence  of  these  lines 
always  allows  a  standardization  to  be  made  as  well  as  an  estimate 
of  the  precision  with  which  the  energy  variable, can  be  established. 

As  an  example  the  data  set  for  Fig.  2b  will  be  used.  Each 
of  ten  experimental  runs  produced  an  independent  value  of  20  for 
each  standard  line.  The  ten  values  of  20  were  then  averaged  to 
obtain  one  value,  listed  in  Table  1  with  the  probable  error  deter¬ 
mined  from  the  replicate  experiments.  Then  from  Bragg's  law  and^ 
the  tabulated  values(13),  2d  (twice  the  interplanar  spacing)  for 
the  mcnochrcmato-  was  calculated  for  each  1  .ne  and  a  final  average 
obtained.  Following  this  procedure  for  every  set  of  data  removes 
the  requirement  for  corrections  involving  the  diffraction  process 
such  as  a  temperature  correction  for  the  monochromator  crystal,  ^ 

J,orentz-polarization  correction,  refraction  and  various  errors 
due  to  misalignment.  The  final  average  for  the  LiF(200)  monochro¬ 
mator,  2d  -  4.0247  ±  3,  is  a  standard  deviation  of  75  ppm  (typical 
of  our  experiments)  which  would  be  ±  0.7  eV  in  the  region  of  the 
Cu  K-edge  of  Fig.  2.  This  value  of  2d  is  used  in  a  computer  program 
to  calculate  (and  plot  a-  in  Fig.  2)  the  value  ot  E  associated 

with  each  data  pair  (IQ,  I)  since 

E  -  12398.04  e\ . 

2d  Sin  0 

Note  that  the  error  reported  for  the  reference  lines  (~40  ppm) 
and  the  precision  to  which  they  can  be  measured  (~50  ppm)  are  comparable 
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and  are  combined  in  the  final  determination  of  2d(  **±  75  ppm). 

There  would  be  little  advantage  in  a  more  accurate  spectrometer 
unless  a  more  accurate  reference  could  be  established. 

m 

Coincidence  Losses  ‘ 

In  Fig.  2a  artifacts  appear  in  the  absorption  curve  at  the 
position  of  the  strong  W  lines  W  L  0^  and  W  t  This  effict  is 

due  to  coincidence  losses  (i  onlinearity)  occurring  in  the  x-riy 
detecting  system.  Although  predictable,  these  features  cause 
serious  errors  in  EXAFS  data  analysis.  An  obvious  solution  is 
to  use  a  very  fast  system  with  a  large  dynamic  range;  however 
instances  will  occur  where  corrections  are  still  necessary. 

We  have  used  the  components  with  the  performance  capabilities 
described  in  Table  2. 

It  is  possible  to  obtain  good  data  with  any  of  these  systems 

(14) 

provided  corrections  are  made.  In  the  method  of  Short  the 

absorption  of  a  single  absorber  is  measured  as  a  function  of  counting 

rate.  When  the  system  becomes  nonlinear  the  measured  absorption 

will  decrease  since  |jx  ■  In  Iq/I  and  Iq,  the  greater  counting  rate, 

» 

will  have  a  proportionally  greater  nonlinearity.  The  performance 
of  the  three  detecting  systems  is  shown  in  Fig.  6a,  b,  c.  Short's 
technique  can  be  used  to  evaluate  a  correction  factor.  This  factor 
is  plotted  in  Fig.  6d.  We  fitted  this  curve  with  an  equation  and 
during  computer  processing  all  data  was  corrected  to  remove  the 
effect  as  in  Fig.  2b  and  c.  The  correction  curve  must  be  evaluated 
to  an  accuracy  comparable  to  the  data,  i.e.  0.1Z. 


I 
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Discussion 

The  experimental  requirements,  of  die  EXAFS  technique  are  within 
range  of  anyone  who  can  dedicate  a  commercial  x-ray  diffractometer, 
power  supply,  and  counting  rack  to  this  use.  The  necessity  for 
multi-pass  averaging  to  obtain  good  statistics  for  the  most  difficult 
(mH  most  interesting)  absorbers  requires  a  full-time  operacion  with 
approximately  24  hours  for  each  pass  through* the  EXAFS  region.  Since 
the  experiment  is  easily  automated  this  is  not  a  handicap.  Indeed, 
our  experiments  are  turned  off  only  for  maintenance  every  few  months. 

The  additional  equipment  items  which  are  essential  are  1)  the 
improved  diffraction  crystal  holder.  (The  intensity  gained  by  opti¬ 
mizing  alignment  through  narrow  slits  can  be  a  factor  of  10.)  -2)  The 
cryostat  and  sample  movement  device.  3)  A  card  punch  for  data  readout. 

Other  methods,  such  as  paper  tape  or  magnetic  taps,  lack  .the 
advantage  of  immediate  visibility  and  easy  correction  of  occasional 
erroneous  data  points. 

After  the  data  is  acquired  and  the  multipass  average  is  obtained 
the  very  considerable  problem  of  data  analysis  remains  which  is  the 
subject  of  the  following  paper. 
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Abstract 

The  technique  of  computer  processing  of  extended  x-ray  absorption 
fine  structure  (EXAFS)  for  the  analysis  of  non- crystalline  solids  is 
discussed.  The  topics  include:  normalization  of  the  data  to  a  per  atom 
scale,  Fourier  filtering  to  remove  the  sloping  background,  and  a  Fourier 
transform  of  the  data  into  a  radial  structure  function.  Particular  factors 
peculiar  to  EX/FS  are  emphasized;  e.g.  upper  and  lower  limits  of  integration, 
the  effect  of  the  phase  shift  on  the  determination  of  interatomic  distances, 
the  peak  shape  and  the  determination  of  coordination  numbers.  The  EXAFS 
spectrum  of  ci/stalline  Ge  is  used  as  an  example. 
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Introduction 


The  previous  paper  described  the  experimental  requirements  for 
obtaining  EXAFS  data.^  This  paper  discusses  the  processing  and  analysis 
of  EXAFS  data  concentrating  on  the  important  factors  which  are  specifi¬ 
cally  related  to  EXAFS.  The  analysis  and  discussion  are  limited  to 
K-shell  absorption  spectra  only.  The  techniques  and  information  involved 
in  L-shell  spectra  will  be  published  elsewhere. ^  The  general  subject  of 
finite  Fourier  transforms  has  been  much  discussed  by  others  involved  in 
x-ray  radial  distribution  analysis.  This  paper  emphasizes  the  processing 
of  the  rw  data,  the  removal  of  the  EXAFS  from  the  average  absorption  back¬ 
ground  using  a  Fourier  filtering  technique,  and  the  normaliration  of  the 
EXAFS  according  to  the  theoretical  expression  so  that  all  data  are  on  the 
same  per  atom  scale.  Fourier  transformation  using  the  point-scattering 
theory  results  in  a  radial  structure  function  centered  on  the  absorbing 
atom.  This  feature  permits  evaluation  of  unique  correlation  functions 
for  each  atomic  species  in  a  multi- component  material  and  is  the  major 
achievement  of  the  EXAFS  technique.  The  particular  problems  of  changing 
from  an  energy  to  '  wave  vector  scale,  selecting  the  upper  and  lower  limits 
of  integration,  the  phase  shift,  and  the  peak  shapes  are  discussed  in  detail. 

Initial  Data  Processing 

The  data  collected  from  our  EXAFS  apparatus  consists  of  the  time.. 

to  collect  a  preset  number  of  counts  (usually  10^)  with  the  absorber  in  and 

out  of  the  beam,  respectively  (T,Tq)  along  with  the  corresponding  step 

number  at  each  angular  position  of  the  spectrometer.  A  typical  data  deck 

consists  of  approximately  200-300  cards  with  4  Pets  of  T,T  and  the 

0 

running  number  per  card.  The  preliminary  data  processing  program  (1)  tab¬ 
ulates  the  initial  data  and  calculates  thw  x-ray  energy,  the  kinetic 
energy  of  the  photoelectron,  and  the  totil  absorption;  (2)  removes  the 
oscillatory  part  of  the  x-ray  absorption  from  the  smooth  monotonic  back¬ 
ground;  and  (3)  normalizes  the  data  to  a  per  atom  scale.  These  three 
program  functions  are  described  below. 
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In  addition  to  tabulating  the  input  data  the  program  routinely 
calculates  the  absorption  and  energy  scales  from  input  data.  Speci- 
fically,  the  total  absorption  is  calculated  at  each  point  using 


yx  -  ln(T/T  ) 
o 


(1) 


as  described  in  the  previous  paper.1  The  x-ray  energy  at  each  point  is 
found  from 


12398.04 

hV  “  idTSTT 


(2) 


where  (2d)  is  calculated  as  described  previously1  and  0  at  each  spectro 
meter  setting  is  found  from  the  initial  angle  and  fixed  angular  step. 
The  kinetic  energy  of  the  ejected  photoelectron  is  then 


E  •  hv-E^  (3) 

where  is  the  initial  binding  energy  of  the  Is  electron  which  was  ejected 
during  the  absorption  of  he  x-ray.  The  tabulated  binding  energies  of 
Bearden  and  Burr  were  used.** 

The  next  step  in  the  processing  involves  removing  the  K  shell 
contribution  for  the  atomic  species  of  interest  from  all  other  absorption 
components  including  other  shells  in  the  same  atom,  absorption  by  other 
atoms  in  the  sample,  and  absorption  by  the  specimen  support  which  is 
usually  A1  foil.  In  the  energy  region  of  interest  these  other  components 
are  usually  far  from  any  characteristic  edges  and  may  be  adequately  des- 
cribed  by  a  Victoreen  formula  of  the  form  y  **  CX  +PX  where  X  is  the  x-ray 
wavelength  and  C  and  D  are  found  by  fitting  this  formula  to  the  measured 
absorption  curve  over  an  energy  range  (typically  a  few  hundred  eV)  Just 
below  the  K-edge  of  interest.  This  fitted  curve  is  then  extrapolated  to 
energies  above  the  absorption  edge  to  approximate  all  of  the  other  absorp¬ 
tion  components  and  is  subtracted  from  the  experimental  curve.  The 
remainder  is  the  K  shell  contribution  to  the  absorption  including  the 
EXAFS  (y„(k)  ). 

k 


2 


The  removal  of  the  oscillatory  part  of  the  x-ray  absorption 
coefficient  (EXAFS)  from  the  smooth,  monotonically-decreasing  background 
absorption  is  a  difficult  problem  since  neither  calculated  nor  empirical 
curves  exist  which  are  accurate  in  the  vicinity  of  the  absorption  edge. 

Rather  than  fitting  some  empirical  function  far  above  the  edge  and  extra¬ 
polating  back  to  the  edge  where  large  uncertainties  may  occur  we  have 
chosen  to  use  a  Fourier  filtering  technique.  This  technique  uses  a  fast 
Fourier  transform  algorithm  on  the  raw  data  (which  is  equispaced  in  0) 
and  takes  2048  data  points  and  analyzes  them  in  terms  of  the  real  and 
imaginary  components  of  1024  frequencies.  The  lowest  frequency  has  a 
wavelength  of  twice  the  range  of  the  data  and  the  remaining  frequencies 
are  the  next  1023  harmonics  of  this  fundamental.  This  spectrum  has  3  com¬ 
ponents:  (1)  contributions  at  low  frequencies  due  to  the  mono  tonically 

decreasing  absorption  coefficient;  (2)  the  EXAFS  signal  at  intermediate 
frequencies;  and  (3)  small  components  from  the  random  statistical  counting 
noise  uniformly  distributed  throughout  the  entire  frequency  range.  The 
lowest  EXAFS  frequencies  were  estimated  on  the  basis  of  the  typical  near 
neighbor  distances  and  were  well  separated  from  the  low  frequency  compon¬ 
ents  which  could  then  be  filtered  by  zeroing  the  first  few  low  frequencies 
and  retransforming  the  spectrum.  This  is  shewn  in  Figure  1  where  70  of 
the  first  512  coefficients  of  a  typical  transform  are  shown.  The  remaining 
coefficients  are  negligible.  The  arrcw  indicates  where  a  typical  cut 
would  be  made  to  filter  out  the  mono  tonic  background.  This  transformed 
spectrum  contains  only  contributions  from  EXAFS  and  the  noise.  The  smooth 
absorption  background  may  be  found  by  subtracting  this  filtered  spectrum 
from  the  raw  data. 

The  entire  normalization  procedure  is  shown  in  Figure  2. 

Figure  2(a)  shows  the  Ge  K-edge  EXAFS  for  amorphous  GeSe  along  with  the 
extrapolated  Victoreen  curve  (dashed  line)  which  was  fitted  to  the  data 
in  the  energy  region  below  the  edge.  In  Figure  2(b)  the  EXAFS  which  was 
obtained  from  the  data  in  Figure  2(a)  using  the  Fourier  filtering  technique 
is  shown  and  in  figure  2(c)  the  smooth  background  curve  obtained  by  sub¬ 
tracting  the  curve  in  Figure  2(b)  from  the  curve  in  Figure  2(a)  is  also 
shown,  figure  2(c)  also  compares  this  smooth  background  curve  V°K0O  to 
both  a  Victoreen  formula7  and  the  empirical  X2*73  curve  suggested  by 
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Figure  1.  Example  of  a  frequency  spectrum  of  the  raw  data  using  a 
fast  Fourier  transform.  Hie  separation  of  EXAFS  from 
the  low  frequency  slope  components  is  apparent.  The 
arrow  indicates  where  the  low  frequency  cutoff  is  made. 
Only  70  of  the  512  components  are  shown;  all  higher 
components  are  negligible. 
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Figure  2a. 


Figure  2b. 


Figure  2c. 


The  measured  K-edge  of  Ge  in  amorphous  GeSe  where  yx 
(absorption  coef.  times  thickness)  is  plotted  ys  the 
photoelectron  kinetic  energy,  E.  The  dashed  line  is  th 
Victoreen  formula  extrapolation  due  to  absorption 
other  than  the  Ge  K-shell. 

x  vs  E  normalized  and  extracted  from  Figure  2a  as  explained 
in  the  text. 

The  smooth  K-shell  absorption  background  vM  vs  E  which 
was  used  to  normalize  x  according  to  equation  (4)  . 
u  (x)  is  obtained  from  Figure  2a  by  subtracting  the 
Victoreen  extrapolation  (dashed  line)  and  Fourier ^  *  8 

to  remove  the  EXAFS.  The  points  6  are  extrapolations  to 
the  edge  from  a  region  1000-1500  eV  above  the  edge  of  8 

both  a8Victoreen  form  and  the  empirical  formula  of  McMaster. 
Both  are  identical  and  approximately  5%  below  the  measured 

curve. 
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McKaster®.  These  curves  were  normalized  to  the  data  in  the  region 
1000-1500  eV  above  the  edge  and  give  essentially  identical  results.  As 
can  be  seen  the  fit  is  good  to  500—1000  eV  but  is  lower  than  the  data 
curve  by  ,v5%  near  the  edge,  which  would  give  noticeable  curvature  to  the 
EXAFS  (Figure  2(b)  if  they  "had  been  used.  We  have  also  checked  the  sensi¬ 
tivity  of  this  technique  to  the  nunber  of  components  which  were  zeroed, 
and  have  found  no  significant  deviation  from  the  smooth  curve  as  long  as 
the  nunber  of  components  zeroed  was  less  than  the  expected  minimum  EXAFS 
frequency  (i.e.  if  no  more  than  the  first  10  frequencies  are  zeroed). 

The  final  step  in  the  initial,  data  processing  involves  normal¬ 
izing  the  data  according  to  the  definition 

X(k)  -  [vK00/v°(k)]-l  (A) 

where  X0O  is  the  EXAFS,  yR(k)  is  the  total  K-shell  contribution  to  the 
x-ray  sbsorption  coefficient  and  u°(k)  is  the  smooth  monotonically  decreasing 
background  of  the  K  shell.  The  EXAFS  defined  in  this  way  normalizes  all 
data  to  the  s^me  per  atom  scale  which  will  «*.llow  absolute  coordination 
numbers  to  be  determined  (as  described  below).  y^(k)  and  y^  CO  are  deter¬ 
mined  as  described  above  and  the  normalized  EXAFS  X00  is  determined  using 
these  functions  and  equation  (A).  The  data  analysis  program  plots  yx  V3.  E 
and  x  vs •  E  and  punches  a  deck  of  x»E  p«irs  which  is  used  as  the  input  for^ 
the  final  data  analysis. 

Fourier  Analysis  of  EXAFS 

Fourier  analysis  of  the  normalized  data  is  based  on  the  point- 

4  5  9 

scattering  theory  of  EXAFS  for  K-shell  absorption 

X(k)  -  1/k  I  -J-E(rj)  gJ(2k)exp(-20jk2)sin[2krJ+2n(k)J  (5) 

Here  g  (2k)  is  the  backward-scattering  amplitude  from  the  surrounding  atoms 
in  the  jth  shell  of  the  excited  p-state  electron  of  wave  number  k  per  unit 
electric  field.  The  direction  of  £  is  radially  out  from  the  excited  atom 
to  the  center  of  an  atom  in  the  jth  atomic  shell ;  E(r^)  is  the  average 
electric  field  induced  at  r^the  position  of  the  jth  shell  of  atoms)  by 
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the  Ionized  atom  at  r  ■  0.  Nj  is  the  total  number  of  atoms  in  the  jth 
shell  and  is  a  measure  of  the  fluctuations  of  the  actual  positions  of 
the  atoms  in  the  jth  shell  about  their  average  position  r  .  The  exponen¬ 
tial  containing  a2  is  a  Debye-Waller-type  factor. 

Since  EXAFS  can  be  viewed  as  a  sum  of  damped  sine  waves  we  can 
formally  invert  the  data  using  Fourier  analysis.  In  most  monatomic 
na*:er^a^8  we  fi-n<i  that  g(2k)  is  independent  of  r^  to  a  good  approximation 
so  that  equation  (5)  can  be  rearranged  and  transformed  with  respect  to 
sin  (2kr+2n(k)  )  to  give 

j^-sin  (2kr+2n(k)  )  (6a) 

exp(-  (r-rj)2/20j)  +  Ar  (6b) 


<Kr)  “ 


/ 

o 


.W. 

8rj°j 


where  equation  (6b)  is  obtained  from  equation  (6a)  using  equation  (5),  and 
$(r)  is  a  radial  structure  function  which  gives  the  location  of  atoms 
around  the  absorbing  atom  in  terms  of  a  series  of  gaussian-shaped  peaks 
whose  amplitudes  are  proportional  to  .  Ar  contains  a  symmetrical 
exponential  term  which  is  not  used  and  other  terms  much  smaller  than  <f>(r). 

While  the  integration  of  equation  (6)  is  illustrative  to  show 
the  kinds  of  information  available  from  EXAFS  it  ignores  several  factors 
which  must  be  accounted  for  in  the  actual  data  analysis  including  the  finite 
upper  and  lower  limits  of  integration,  the  fact  that  g(2k)  and  n  (k)  are  not, 
in  general,  known  and  that  the  definition  of  k  is  uncertain.  Specifically 
thi  integral  that  is  actually  calculated  is 

k 

4(f)  -  VTT  /  max  x(k)  exp  (2ikr)dk  (7) 

2  tt  k  . 
min 

The  problems  involved  in  equation  (7)  fall  in  2  categories 
(evaluation  and  interpretation)  and  will  be  discussed  separately. 
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A)  Evaluation  -  The  photoelectron  wavevector  is  evaluated  from  the 
definition 


-  E  -  (Uq-V0  (8) 

2m 

where  E  is  defined  by  equation  (3)  and  is  referenced  to  the  Fermi  energy, 

U  is  the  average  inner  potential  measured  from  the  vacuum  level  and  W  is 
o 

the  work  function.  For  metals  Uq-W  is  on  the  order  of  the  Fermi  energy 
measured  from  the  bottom  of  the  conduction  band.  In  semiconductors  and 
insulators  the  zero  of  kinetic  energy  is  not  as  easy  to  define  but  may  be 
approximated  by  averaging  the  pseudo-potential  over  the  energy  range  of 
interest.  In  either  case  errors,  AE,  of  a  few  eV  may  arise  in  the  choice 
of  U  -W  used.  These  errors  do  not  significantly  affect  the  data  since 
even  at  the  energy  corresponding  to  the  lower  limit  of  integration  (typi¬ 
cally  E.  'v-SOeV),  AE/E  is  negligible.  We  have  verified  this  on  several 
spectra  by  varying  Uo-W  over  a  suitable  range  and  found  no  significant 
effects  on  the  transforms. 

Another  problem  in  evaluation  is  the  choice  of  upper  and  lcwer 
limits  of  integration  in  equation  (7).  The  upper  limit  is  determined 
primarily  by  the  experimental  range  of  energies  which  have  been  measured, 
typically,  1000-1500  eV  above  the  edge  (which  corresponds  to  a  k  range  up  ^ 
to  15-20  A  )  unless  other  absorption  edges  appear  in  the  data.  It  is 
usually  desirable  to  have  k  as  large  as  possible  to  reduce  the  termination 
error  found  in  Fourier  transforms  which  are  truncated^  at  finite  k,  although 
these  effects  are  not  the  primary  cause  of  deviations  from  the  gaussian- 
peak  shape  predicted  by  equation  (6b)  as  will  be  discussed  below.  The 
lower  limit  of  integration  is  generally  determined  by  the  requirement  that 
it  is  far  enough  above  the  edge  so  that  the  assumptions  of  the  point 
scattering  theory  are  valid  and  any  band  structure  effects  on  the  density 
of  states  have  disappeared.  There  is  no  precise  point  at  which  this 
occurs  but  is  typically  ^30-50  eV  and  may  be  verified  by  choosing  several 
lower  limits  and  selecting  the  one  which  produces  the  cleanest  transform 
(i.e.,  lew  noise  at  r<l£). 
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The  actual  transform  may  be  done  either  by  direct  integration 
or  by  using  a  Fast-Fourier  algorithm.  In  order  to  process  large  amounts 
of  data  with  the  capability  to  make  real  time  decisions  we  have  developed 
a  program  using  a  display  scope  with  interactive  capability10.  The  out¬ 
puts  are  the  real  (cosine),  imaginary  (sine),  and  magnitude  of  the  trans¬ 
form.  In  Figure  3  the  EXAFS  for  crystalline  Ge  over  the  range  to  be  trans 
formed  is  shown.  In  Figure  4  the  transform  of  the  data  from  Figure  3  is 
shown  over  the  range  of  k  from  4.1  to  15.8a"1,  where  4(a)  is  the  real  part 
4(b)  the  imaginary  part,  and  4(c)  is  the  magnitude  of  the  transform.  As 
can  be  seen  the  imaginary  and  magnitude  have  several  peaks  which  contain 
structural  information.  The  factors  which  affect  the  interpretation  of 
these  transforms  and  the  informati->n  contained  in  the  transform  is 
discussed  in  the  next  section. 

B)  Interpretation  -  There  are  several  qualitative  features  which  help 
understand  the  data.  First,  as  can  be  Jeen  in  Figure  4,  there  is  very 
little  amplitude  in  the  region  below  1A.  This  condition  can  usually  be 
achieved  by  proper  selection  of  the  lower  limit  of  integration  and  indi¬ 
cates  that  there  are  no  large  normalization  or  background  subtraction 
errors  which  have  been  known  to  affect  the  transforms  strongly  in  this 
region  ’  .  Because  of  the  uncertainties  which  have  beer,  discussed  in 

renoving  the  EXAFS  from  the  background  absorption  it  is  likely  that  some 
error  has  been  introduced  in  the  transform;  however,  our  experience  is 
that  these  errors  can  be  made  small.  - 

The  relationship  between  the  real  part,  imaginary  part  and 
magnitude  of  the  transform  is  apparent  in  the  first  large  peak  in  the 
data  (between  1-3A).  If  the  theoretical  expression  from  equation  5  is 
inserted  into  equation  7  and  it  is  assumed  that  n=0  and  the  k  dependence 
of  the  factor  g  (2k)/k  is  neglected  the  result  would  be  a  symmetric 
gaussian-shaped  imaginary  part  and  an  antisymmetric  real  part  centered 
on  the  shell  distance,  r  ^ .  This  is  what  is  observed  in  Ge  for  the  first 
shell.  The  peak  shape  of  the  imaginary  part  has  deep  negative  side 
lobes,  which  will  be  discussed  later,  but  this  does  not  change  the  basic 
relationship  between  the  two  parts  of  the  transform.  In  fact,  the  two 
parts  of  the  transform  contain  complementary  information  so  that  only 
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Figure  4a, 
4b,  4c. 


The  real  part,  imaginary  part,  and  magnitude,  respectively, 
of  the  Fourier  transform  of  the  EXAFS  in  Figure  3  plotted 
$(r)  vs  r(X) .  The  numbers  over  the  peaks  locate  the  first 
6  coordination  shells.  The  <t>(r)  scale  in  4c  is  twice  that 
of  4a  and  4b. 
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one  part  need  be  studied.  The  situation  is  eomewhat  complicated  by 
the  finite  constant  phase  shift.  By  using  simple  trigonometric  iden- 
lltlti  it  can  be  shown  that  the  real  and  Imaginary  parts  now  are 
admixtures  of  the  q*0  transforms  so  that  neither  part  is  completely 
symmetric  or  antisymmetric  unless  the  phase  is  a  multiple  of  tt/2. 
Location  of  the  peak  position  may  still  be  made  using  the  magnitude  of 
the  transform  which  is  independent  of  any  constant  phase  shift  and  will 
be  peaked  at  r^r^.  In  the  case  of  Ge  it  appears  that  the  constant  part 
of  the  phase  is  approximately  zero  since  the  real  and  imaginary  parts 
of  the  Fourier  transform  have  the  shapes  described  above.  Note  that 
the  magnitude  (Figure  Ac)  is  peaked  at  the  same  place  as  the  imaginary 
(Figure  A—)  part  but  is  smoother  and  broader  which  in  some  ways  makes 
it  easier  to  identify  shell  locations  since  it  reduces  the  effects  of 
the  side  lobes. 

To  locate  the  peak  position  accurately  it  is  necessary  to 
understand  the  phase  shift  more  fully.  Estimates  of  the  phase  shifts 
using  the  method  have  shown  that  for  most  elements  the  phase 

shift  is  approximately  linear  in  k  over  the  EXAFS  energy  range  with  a 
negative  slope,  i.e. 


rj(k)  -  -ak+6  (9) 

-  r 

with  a,  8  constant.  This  linearity  has  been  verified  in  the  case  of  Cu 

using  a  more  exact  calculation  of  phase  shifts  using  a  Hartree-Fock 

potential.^  Substituting  equation  (9)  into  equation  (5)  gives  the 

argument  for  the  sine  wave  of  (2kr-2ak+2B) .  The  effect  of  the  constant 

phase  26  has  been  discussed  above;  however,  as  first  pointed  out  by  Mott1^, 

the  linear  k  dependent  term  shifts  the  "frequency"  of  the  sine  wave  from 

Tj  to  rj-a.  In  the  Fourier  transform  this  has  the  effect  of  shifting 

all  peaks  toward  the  origin  by  an  amount  a.  In  crystalline  materials 

where  the  distances,  r^,  are  known  the  transforms  may  be  used  to  evaluate 

a.  This  is  shewn  in  Table  I  where  the  data  for  the  transforms  shown  in 

Figure  A  is  tabulated.  Column  five  of  the  table  shews  that  for  the  first 

o 

3  shells  where  the  data  is  most  reliable  a  value  of  a  -  .28A  may  be  used. 
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Table  I.  EXAFS  Structural  Data  For  Ge 


j 

0 

r^(A) 

»  o 
ri(A) 

0 

a  (A) 

y*i 

1 

4 

2.45 

2.17 

.28 

l 

2 

12 

4.00 

3.72 

.28 

.45 

3 

12 

4.68 

4.41 

.27 

.30 

4 

6 

5.65 

5.38 

.27 

.04 

5 

12 

6.16 

5.83 

.33 

.06 

6 

24 

6.92 

6.60 

.32 

.13 

r^  ■  known  jth  shell  distance  for  dia.  cub.  Ge  (d=5.65A) 
r  -  measured  peak  position  from  Figure  4b 

J 

°  ‘  Yrj 

» 

Aj/A^  *  amplitude  of  the  imaginary  transform  at  r^  normalized 
to  the  first  shell  amplitude 
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This  data  also  shows  at  least  6  shells  in  Ge  may  be  located  although 

o  o 

there  are  other  features  in  the  data  at  ^5A  and  ^6.2A  which  can  t  be 
explained  by  the  simple  point  scattering  model.  The  value  of  a  which 

o 

can  be  determined  from  kncwn  structures  is  typically  accurate  to  +.02A 
which  gives  a  total  accuracy  of  2-3%  in  the  distance  measurements.  To 
determine  distances  in  materials  where  the  crystal  structure  is  not 
known  (such  as  amorphous  materials  or  biological  systems)  assumptions 
must  be  made  about  the  phase  shift.  Generally,  the  phase  shift  is  domi¬ 
nated  by  the  potential  well  created  by  the  inner  vacancy  in  absorbing 
atom  itself  and  not  by  its  surroundings  so  that  all  Ge  atoms  should  have 
approximately  the  same  phase  shift  regardless  of  the  compound  it  is  in. 

This  is  verified  by  the  data  in  Table  II  where  the  k  dependence  of  the 
phase  shifts  in  a  variety  of  materials  whose  structures  are  known  are 
listed.  The  Ge  phase  shifts  in  the  Ge  compounds  agree  very  well  with  the 
value  of  a  -  .28  determined  from  crystalline  Ge.  The  values  of  a  for 
several  other  elements  in  the  third  row  of  the  periodic  table  are  listed 
and  show  that  a  is  between  .25  and  .41  in  this  region  of  the  periodic 
table.  To  minimize  the  uncertainties  in  assigning  a  for  unknown  structures 
we  have  used  standard  compounds  whose  composition  and  local  structure  are 
known  and  whose  local  environment  is  expected  to  be  similar  to  the  local 
structure  being  investigated.  As  a  general  rule  we  have  found  that  com¬ 
pound?  among  atoms  in  the  same  row  of  the  periodic  table  (especially  if  ^ 
the  valence  of  the  absorbing  atom  is  small)  tend  to  have  phase  shifts  similar 
to  the  pure  element  whose  edge  is  being  studied.  We  have  also  found  that 
in  some  systems,  particularly  certain  oxides  with  ionic  bond,  there  are 
significant  deviations  from  the  elemental  phase  shifts.  Further  experi¬ 
mental  and  theoretical  studies  are  underway  in  order  to  quantitatively 
understand  the  observed  phase  shifts  and  their  origins  and  make  reliable 
structural  measurements  in  a  wider  variety  of  systems. 

As  mentioned  above  the  shapes  of  the  peak  differ  significantly 
from  the  gaussian  shape  predicted  by  equation  (6b)  even  if  termination 
effects  are  taken  into  account.  This  is  because  in  the  actual  transform 
(equation  (7))  the  k  dependence  of  the  factor  g(2k)/k  has  not  been  removed 

from  the  data  so  that  the  experimental  transform  is  a  convolution  of  the 
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Table  II.  Phase  Shift  k  Dependence  (a  from  equation  9) 
for  Several  Elements  and  Compounds 


Mi  terial 

a 

Material 

a 

Cr 

.36 

Ge 

.28 

Fe 

.26 

Se 

.26 

Ni 

.25 

(Ge)02 

.27 

Cu 

.38 

(Ge)Se2 

.29 

Zn 

.41 

(Ge)Se 

.29 

As 

.27 
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result  of  equation  (6b)  with  the  Fourier  transform  of  g(2k)/k.  The 
general  shape  of  this  function  is  the  envelope  of  the  EXAFS  shown  in 

o-l 

Figure  3  which,  for  many  systems,  increases  to  a  peak  about  k  ■  6-7  A 

and  decreases  at  higher  k.  This  envelope  also  includes  temperature 

effects  but  these  decrease  monotonically  with  k  and  do  not  give  the  peaking 

which  is  observed.  The  transforms  of  a  variety  of  similar  envelope 

functions  have  been  studied  by  Waser  and  Schomaker  7  who  show  that  the 

results  are  transforms  with  large  negative  side  lobes  as  observed  in 

Figure  A.  At  present  the  functional  form  of  g(2k)  is  not  known  well 

enough  to  separate  its  effect  from  the  data.  This  complicates  the 

analysis  particularly  at  large  r  where  the  shells  are  close  enough 

together  so  that  the  side  lobes  interfere  with  adjacent  peaks  making 

peak  identification  uncertain.  The  peak  shape  also  makes  interpretation 

of  the  amplitude  difficult.  From  equation  (6b)  the  amplitude  of  each 

2 

peak  is  proportional  to  N,E(r  )/o,r.,  and  by  a  standard  integration  the 

j  j  j  j  2 

area  under  the  peak  is  proportional  to  NjE(rj)/r^.  However,  the  area  is 
difficult  to  measure  accurately  so  the  amplitude  is  usually  used  to 
cbtain  any  information  about  number  of  atoms  per  shell.  Since  r^  is  deter¬ 
mined  from  the  peak  location,  the  amplitude  contains  information  propor¬ 
tions]  to  NjE(rj)/Oj.  When  the  data  is  normalized  according  to  equation  A, 
all  data  are  on  the  same  scale  and  could  be  compared  directly  to  deter¬ 
mine  the  coordination  numbers,  N^,  if  E(r^)/aj  were  known.  Normally,  this, 
is  not  the  case  so  reliable  coordination  numbers  are  obtained  by  empirical 
comparisons.  In  some  crystalline  systems  may  be  calculated  and  the 
measured  amplitude  used  with  the  known  crystal  structure  to  determine 
E(rj  );  however,  it  is  usually  not  possible  to  measure  enough  amplitudes 
reliably  to  differentiate  between  possible  forms  for  E(r^)  (e.g.  power  1®? 
vs.  exponentail) .  It  is  also  possible  to  measure  the  temperature  dependence 

of  the  amplitudes  to  determine  c.  experimentally.  In  multi-component 

J  17 

materials  where  only  one  absorption  edge  can  be  measured,  we  have  assumed 
that  E(r^)  is  primarily  a  function  of  local  environment  and  used  standard 
compounds  whose  environment  and  coordination  are  expected  to  be  similar  to 
the  unknown  system  to  determine  N^.  We  estimate  that  this  procedure  can 
give  coordination  numbers  which  are  good  to  ^20%  for  the  first  few  shells. 
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Recently,  we  have  extended  this  procedure  to  making  absolute  measure¬ 
ments  of  the  first  neighbor  coordination  of  both  species  in  binary  com- 
pounds  \diere  each  atom's  EXAFS  can  be  measured.  It  was  possible  to 
determine  the  number  of  each  kind  of  first  neighbor  atom  around  both 
species.  This  technique  was  limited  to  systems  for  which  good  standards 
could  be  used  for  a  reference,  but  did  give  reliable  results  for  amorphous 
GeSe2  and  As^Se^  which  are  similar  to  their  crystalline  polymorphs.  To 
obtain  more  reliable  coordination  numbers  will  require  a  better  under¬ 
standing  of  the  meachanisms  which  determine  E(r)  and  g(2k),  so  that  these 
factors  can  be  determined  reliably  without  the  use  of  standard  compounds. 
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Pis  cuss  ion 


The  Fourier  transform  of  EXAFS  data  can  be  used  to  determine 
the  local  atomic  environment  of  each  kind  of  atom  in  multi-component 
non-crystalline  materials.  It  is  a  general  and  powerful  tool  for  study¬ 
ing  any  physical  state  of  matter,  solid,  liquid,  or  gas.  It  has  been 

5  9  18  17 

applied  to  amorphous  semiconductors,  ’  ’  supported  catalysts  and 
biological  materials  (yet  to  be  published). 

The  data  processing  is  concerned  with  extraction  of  the  EXAFS 

from  data  containing  the  effects  of  other  physical  processes  in  a  form 

suitable  for  Fourier  analysis,  which  is  a  difficulty  common  to  other 

scattering  methods.  However,  there  are  particular  factors  peculiar  to 

EXAFS  such  as  the  phase  shift,  definition  of  the  wave  vector  scale,  and 

an  unknown  scattering  factor  which  are  the  primary  concerns  of  this 

paper.  We  have  addressed  these  problems  in  a  completely  pragmatic  way 

by  using  well-matched  standard  materials  to  determine  these  parameters 

and  then  applying  them  to  the  unknown  materials.  As  such,  the  technique 

is  useful  for  studying  the  atomic  arrangement  in  non- crystalline 

materials  but  suffers  a  lack  of  completeness.  As  more  materials  are 

studied  and  the  parameters  become  more  accurately  defined,  we  trust  that 

theory  will  catch  up.  In  addition  to  the  structural  information 

emphasized  here,  the  x-ray  absorption  edges  contain  a  wealth  of  physics 

which  is  becoming  apparent  as  we  have  tried  to  understand  the  EXAFS 

technique;  e.g.  some  L-spectra  show  giant  resonances  at  the  edge  with 

2 

a-d  mixing  in  the  EXAFS  and  the  role  of  impurity- charge  shielding  in 

9 

EXAFS  has  been  discussed. 
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ABSTRACT 


A  computer  program  Is  described  which  allows  real-time 
data  analysis  decisions  by  means  of  a  display  scope  for 
Fourier  transforms,  data  manipulation,  and  Inverse  trans¬ 
forms  of  x-ray  data.  The  example  given  analyzes  the  Ge 
K-edge  extended  x-ray  absorption  fine  structure-,  however 
anv  x-ray,  electron,  or  neutron  scattering  data  could  be 


INTRODUCTION 


Fourier  transforms  of  experimental  data  are  widely  used  in  many  areas 
of  physics  including  x-ray  and  electron  scattering  and  extended  x-ray 
absorption  fine  structure  (EXAFS)  experiment:  as  well  as  many  other 
applications.  Often  these  transforms  involve  large  amounts  of  data 
and  are  very  time-consuming;  however,  the  development  of  the  Fast 
Fourier  Transform  (FFT)  algorithm  (Cooley  &  Tukey,  1965)  resulted  in 
substantial  time  savings.  Data  analysis  may  still  be  time-consuming, 
particularly  if  real-time  decisions  must  be  made  regarding  optimum 
termination  points,  multiplication  by  termination  functions,  inverse 
transforms,  and  display  of  final  data,  which  traditionally  means  several 
individual  computer  runs  with  resultant  turnaround  time. 

These  shortcomings  have  been  overcome  with  the  development  of  a  data 
analysis  program  which  uses  a  vector  display  scope  on  line  with  a  large 
digital  computer.  This  interactive  computer  program  allows  rapid  initial 
processing,  real-time  decisions  regarding  the  data  with  a  light. pen  or 
typewriter,  and  FFTs  all  of  which  require  only  a  few  minutes  per  data  case. 
In  the  following  sections  the  hardware,  software  and  a  sample  data  analysis 
of  EXAFS  will  be  used  to  demonstrate  the  program's  capabilities 

The  program  has  been  developed  to  analyze  digital  EXAFS  data  (see  Sayers, 
Lytle,  &  Stern,  1971  for  examples),  but  can  be  used  for  any  other  kind  of 
Fourier  analysis  problem.  Data  in  graphical  form  may  also  be  analyzed  by 
digitizing  first.  Most  scattering  experiments  (x-ray,  electron,  neutron 
or  EXAFS)  are  conducted  so  that  data  are  accumulated  in  equi -spaced  angular 
increments  which  are  not  equally  spaced  in  the  transform  variable.  This 
interpolation  is  accomplished  as  a  first  step  in  this  program,  hence  any 
data  equi -spaced  or  not  may  be  analyzed. 


44< 


1.  COMPUTER  HARDWARE  * 


The  hardware  for  this  program  is  a  System/360  model  65  IBM  computer, 
attached  directly  to  an  IBM  2250  vector  display  scope  with  a  32-key 
function  keyboard,  a  standard  typewriter  type  entry  keyboard,  and  a 
light  detecting  pen.  The  scope  is  similar  in  appearance  to  a  TV 
screen.  The  light  detecting  pen  can  be  touched  directly  on  the 
screen  yielding  the  precise  coordinates  of  contact,  which  can  be 
used  both  for  terminating  the  data  and  issuing  program  instructions. 
Internal  core  storage  is  limited  to  280K  bytes  at  our  facility  because 
of  the  multi -processing  requirements  of  the  computer.  This  problem 
has  been  reduced  by  development  of  an  FFT  algorithm  (Singleton,  1967) 
which  uses  a  relatively  small  amount  of  in-core  storage  and  uses  10 
auxiliary  mass  storage  units  for  data  storage.  The  input  data  is 
entered  on  cards  with  the  program  residing  on  permanent  disc  storaye. 
Printed  output  goes  on  a  standard,  high-speed  printer.  Graphical 
output  (hard  copies)  of  any  screen  display  may  be  obtained  on  either 
SC4020  or  CALCOMP  plotting  hardware  at  the  user's  option.  This  output 
Is  written  on  two  7-track  tapes  to  be  processed  offline.  The  entire 
program  requires  280K  bytes  of  internal  mass  storage  and  operates  well 
In  the  multi-variable  task  processing  environment. 

2.  COMPUTER  SOFTWARE 

The  Interactive  computer  graphics  program  consists  of  one  long  driver 
program  that  calls  numerous  subroutines  from  the  IBM  2250  graphics 
subroutine  package  (IBM,  1968)  and  the  FFT  subroutine.  All  mathematical 
manipulation  with  the  exception  of  the  transform  is  performed  by  the 
main  program.  The  IBM-supplied  software  facilitates  communication  with 
the  main  program  in  controlling  light-pen  and  function  key  interrupts, 
display  regeneration,  and  typed  entries. 

At  our  facility  core  storage  for  scope  jobs  is  limited  to  280K  bytes  of 

13 

core  which  represents  a  problem  since  Fourier  transforms  of  2  complex 
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2.  (Continued) 


data  pairs  must  be  performed.  The  graphics  subroutine  package  takes 
about  170K  bytes  of  storage,  leaving  only  110K  bytes  for  the  main 
program  and  the  FFT  algorithm.  Generally,  FFT  algorithms  (Bingham 
&  Godfrey,  1967)  need  at  least  2.5  times  the  size  of  the  transform 
for  internal  storage  in  order  to  perform  an  in-core  transform.  Since 
this  amount  of  storage  was  simply  not  available  for  processing  in  a 
graphics  mode,  an  algorithm  was  chosen  that  uses  auxiliary  mass  storage 
(Disc  Files)  (Singleton,  1967).  This  FFT  algorithm  occupies  only  16K 
bytes  of  storage  plus  20K  bytes  for  high-speed  buffers  to  perform  a 
double  precision  transform.  All  of  the  routines  used  including  the 
FFT  subroutine  were  coded  in  FORTRAN  with  the  exception  of. certain 
GSP  routines. 

3.  ANALYSIS  PROCEDURE 

In  order  to  demonstrate  the  capabilities  of  this  program,  a  sample  data 
examination  is  presented  using  the  K-euge  EXAFS  from  crystalline 
germanium  (Ge).  EXAFS  is  the  series  of  oscillations  appearing  on  the 
high  energy  side  of  x-ray  absorption  edges.  The  theory  of  EXAFS  (Sayers, 
et  al.,  1971)  is  given  by  •  ^ 

x(k)=I  NJ  g(k)  exp(-c^k2)  S1n(2krj  +  2n) 
j=l  7? 

where  the  sum  is  over  the  shells  of  Nj  atoms  surrounding  the  absorbing 

atom  at  a  distance  r.,  g(k)  is  the  energy  dependence  of  the  scattering 

of  the  photo-electron,  the  exponential  is  a  Debye-Waller  type  temperature 
term  and  n  Is  a  phase  shift  introduced  by  the  absorbing  atom.  Because 
of  the  basic  sine  wave  dependence  of  x(k)  we  have  shown  that  Fourier 
analysis  of  the  data  gives  peaks  located  at  the  atomic  distances  whose 
areas  are  proportional  to  Nj.  This  technique  has  bden  applied  to 
studies  of  amorphous  semiconductors  (Sayers,  et  al . v  1971),  catalysts 
{Lytle,  et  al.,  1974),  and  charge  shielding  in  metals  (Stern  S  Sayers, 
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1.  DISPLAY  CHI  VS  K  FOR  THE  DATA  OUST  READ  IN. 

2.  CHOOSE  LIMITS  IN  K  SPACE  FOR  FOURIER  TRANSFORM. 

3.  DISPLAY  CHIBAR  VS  KBAR  FOR  THE  RANGE  TO  BE  TRANSFORMED 
AFTER  THE  DATA  HAS  BEEN  EQU I SPACED  IN  K. 

4.  MULTIPLY  BY  CONVERGENCE  FACTOR,  A*EXP(-B* k2) . 

5.  PERFORM  THE  FOURIER  TRANSFORM. 

6.  DISPLAY  THE  TRANSFORM: 

REAL 

IMAGINARY 

MAGNITUDE 

7.  MULTIPLY  THE  TRANSFORM  BY  CrN. 

8.  SELECT  THE  LIMITS  FOR  THE  INVERSE  TRANSFORM. 

•  r  - 

9.  PERFORM  THE  INVERSE  TRANSFORM. 

10.  DISPLAY  THE  INVERSE  TRANSFORM. 

•  • 

11.  READ  THE  NEXT  DATA  CASE. 


FIGURE  1  -  "Menu"  of  Program  Options 
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3.  (Continued) 


One  or  more  decks  of  Chi  versus  E  (E  is  the  energy  above  the  x-ray 
absorption  edge)  are  read  into  the  computer.  The  program  displays 
the  "menu"  (shown  in  Figure  1)  on  the  scope  listing  the  analysis 
options  available.  Normally  analysis  proceeds  step-by-step  return¬ 
ing  to  the  menu  after  each  step;  however,  if  desired,  any  step  may 
be  repeated.  The  options  are  shown  in  the  flow  chart  in  Appendix  I 
and  are  selected  by  touching  the  light-pen  on  the  line  of  the  menu 
desired. 

Normally,  the  first  operation  is  to  display  the  input  data  with  a 
change  of  independent  variable  from  energy  (E)  to  the  electron  wave 
vector  k  which  is  the  appropriate  Fourier  variable  and  is  related 
to  E  by 

1/2 

k  *[fjr  (E-»ojl 

where  m  is  the  electron  mass,  A  is. Planck's  constant,  and  V0  is 
an  inner  potential  which  defines  the  zero  of  kinetic  energy. 

Pressing  either  a  plot  button  (for  CALCOMP  or  SC4020  plots)  or  the 
CONTINUE  button,  returns  the  menu  display.  At  this  point  select 
line  two  of  the  menu  and  type  the  upper  and  lower  transform  limits 
in  k  with  the  standard  keyboard  or  press  continue  and  CHI  versus 
k  will  again  appear  to  allow  selection  of  the  data  range  to  be  trans¬ 
formed  (using  the  light-pen)  or  choose  any  combination  of  typed  limits 
or  light-pen  are  automatically  marked  with  a  "-f-".  This  display  may 
be  plotted  and  again  the  menu  returns  after  either  plotting  or  a 
continue  is  pressed.  The  initial  data  CHI  vs.  k  for  crystalline  Ge, 
with  a  typical  transform  range  selected  by  the  light-pen  is  shown  in 
Figure  2.  ^ 
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3.  (Continued) 

The  FFT  algorithm  requires  input  data  which  is  e^ji-spaced  in  k. 

This  is  generated  by  touching  the  next  line  of  the  menu  with  the 
light-pen  which  uses  a  parabolic  interpolation  routine  to  equi -space 
the  data  (Chibar  vs.  Kbar).  The  number  of  points  displayed  in  the 
range  of  the  transform  is  proportional  to  the  number  of  input  data 
points.  These  points  are  followed  with  zeros  to  provide  a  total  of 
213  complex  data  pairs.  There  are  400-1000  values  of  Chibar  in  a 
typical  transform  range  and  the  remaining  7000  or  so  points  are 
zeros.  This  technique  (Bingham,  et  al . ,  1967)  is  widely  used  to 
achieve  higher  resolution  in  the  transform.  By  extending  the  length 
of  the  period  in  k-space  through  adding  zeros,  the  distance  between 
harmonics  (r-space)  is  reduced.  The  transform  input  files  are  filled 
and  the  display  for  Chibar  vs.  Kbar  (the  equi-spaced  values  of  k), 
appears  on  the  screen.'  The  user  has  an  option  of  multiplying  Chibar 
by  a  convergence  factor  A-exp(-B-k2)  where  default  values  are  A=1 
and  B=0.  A  plot  of  CHIBAR  vs.. KBAR  with  a  convergence  factor 
(02=b=.91 )  is  shown  in  Figure  3.  This  data  may  also  be  printed  on  the 

line  printer. 

The  next  line  performs  the  actual  transform.  While  this  operation  is 
progressing,  the  function  buttons  and  scope  display  are  off  to  save, 
core  storage.  Typical  transforms  take  about  3  minutes  depending  on 
other  multiprocessing  demand.  We  are  currently  improving  this  version 
to  reduce  this  time  to  less  than  1  minute.  When  the  transfonn  is 
complete  the  menu  reappears  on  the  scope  and  the  DISPLAY  TRANSFORM 
option  is  selected.  An  option  is  presented  to  multiply  the  transform 

byCrN  (default  values  A=1,N=0).  The  real  part,  imaginary  part,  and 

magnitude  (magnitude  =  [(real)2  ♦  (imag)2]1/2  of  the  Fourier  transform 
between  r=0  and  10  A  are  consecutively  displayed  by  pressing  the  plot 
or  continue  buttons.  This  data  may  al sc  be  printed.  The  real  and 

imaginary  part  as  well  as  the  magnitude  of  the  transform  of  the  data 
from  Figure  3  are  shown  in  Figure  4  where  Cr  =  ^r2: 


\ 


\ 
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Radial  Distance,  r  (A) 

Figures  4a,  b  and  c  The  real  part,  imaginary  part 

and  magnitude  of  the  FFT  for  the  data  of  Fig.  2.  The 

vertical  scale  is  proportional  to  the  number  of  atoms 

on  each  coordination  shell  and  has  been  multiplied  by 

■r^.  The  horizontal  scale  is  the  radial  distance  from 
o 

the  absorbing  atom  in  A. 


Electron  Wave  Vector,  K(A  I 

Figure  5  The  inverse  FFT  for  the  first  peak  in 
Fig,  4c  for  1.49  <r  <2.81  X  plotted  EXAFS  ampli 
tude  vs  radial  distance,  r. 
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After  completion  of  this  step  the  menu  Is  displayed  again  on  the 
screen.  The  options  on  lines  8-10  allow  all  or  part  of  the  trans¬ 
form  In  r-space  to  be  Fourier  transformed  back  to  k-space.  This  is 
very  useful  in  structural  analysis  and  other  applications  because  it 
allows  the  contribution  of  a  particular  feature  or  range  in  r-space 
to  be  studied  in  k-space  without  interference  from  other  features. 

From  line  8  of  the  menu,  the  range  of  data  to  be  inverse  transformed 
is  selected  using  the  typewriter  or  light-pen  as  was  done  with  the 
initial  transform.  All  data  below  the  lower  limit  in  r  and  above 
the  upper  limit  in  r  are  set  to  zero.  No  interpolation  is  needed 
because  the  data  from  the  original  transform  is  equi-spaced  in  r. 
Selection  of  the  next  line  performs  the  inverse  transform  just  as 
was  done  previously,  and  the  next  line  of  the  menu  displays  the  re¬ 
sult.  Figure  5  shows  the  result  of  isolating  a  peak  in  r-space  (the 
first  peak  in  Figure  4C)  and  subsequent  inverse  Fourier  transformation. 
This  shows  what  part  of  the  signal  Chi  vs.  k  contributed  to  the  selected 
peak  In  r-space.  A  comparison  of  Figure  5  and  Figure  2  shows  which 
features  in  the  original  data  can  be  identified  with  the  first  peak. 

If  no  other  inverse  transforms  are  required,  the  user  has  the  option 
of  taking  another  transform  over  a  different  range  of  the  data,  multiply- 
Ing  by  a  convergence  factor,  or  multiplying  the  transform  by  Cr  or 
going  on  to  the  next  data  case.  The  sample  run  which  has  been  illustr¬ 
ated  here  would  take  less  than  10  minutes  in  actual  time  on  the  computer 
so  that  several  cases  per  hour  can  be  analyzed. 


We  have  found  this  program  invaluable  as  an  aid  for  processing  EXAFS 
data.  This  interactive  program  package  is  available  for  general  use 
at  Boeing  Computer  Services,  Inc.,  Seattle,  Washington.  From  our 
experience  we  have  found  certain  operations  which  are  well  known  in 
advance  and  repetitious  in  nature.  These  problems  ar.p  handled  with  a 
similar  batch  program  which  does  not  require  user  interaction.  This 
program  is  also  aval ’’able. 


*  Research  supported  in  part  by  the  Advanced  Research  Projects  A9ency 

and  monitored  by  the  U.S.  Army  Research  office  under  Contract  DAHC04-71 -C001 0 . 
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